function loglike = estimateLoglikelihood(mix, act)

valid_inds = find(mix.priors > 0);  % find all the valid components
prob = act(:, valid_inds)*(mix.priors(valid_inds))';
loglike = sum(log(prob+realmin));
